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Spectral Flow and Quantum Theory of Dissipation 
in the Vortex Core of BCS superconductors 

Masahiko HayashiB 

Graduate School of Information Sciences, Tohoku University, Sendai 980-8579 
(Received February 1, 2008) 

The dissipation process in two dimensional BCS superconductor due to the quasiparticles in 
the core of moving vortex is studied from a quantum mechanical point of view, especially paying 
attention to the spectral flow of the bound state energies in the vortex core. In order to clarify the 
nature of the spectral flow, we performed a numerical study of a finite system and, by extending 
the analysis to include the effect of impurities, we discuss the quantum mechanical origin of the 
dissipation in the vortex core. 
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§1. Introduction 

Recently the vortices in type-II superconductors in the 
quantum limit, where the level spacing of the Caroli- 
de Gennes-Matricon's (CdGM's) bound states in the 
vortex core, hu>o = A^/Ep (An and Ep axe bulk en- 
ergy gap and Fermi energy, respectively) ,cP is larger 
than the energy scale of thermal fluctuation, ksT (ks 
and T are Boltzmann constant and temperature, respec- 
tively), have been attracting much attention stimulated 
by the study of high T c oxide and organic superconduc- 
tors. Especially in the moving vortices, the effects of 
bound states are believed to be essential in understand- 
ing the vortex dynamics. In order to treat this problem, 
an interesting idea has been proposed, that is the spec- 
tral flow flu If This is based on the finding that the 
bound state energies in the vortex core change continu- 
ously when the vortex is moving. Employing this idea, 
the nondissipative part of the vortex dynamics has been 
studiedJntcnsively, although there still remains contro- 
versy .Em' 

On the other hands, it is also an interesting prob- 
lem to see how the dissipative part of vortex dynam- 
ics is modified from the Bardeen-Stephen's classical the- 
oryEP in the quantum limit. Several attempts have been 
made to clarify this point. Until now, however, most 
of the studies were focused on the transition between 
CdGM'^J^piuLstates due to the interaction with impu- 
ritiesflBli3'EJIiiJ In this letter, we argue that the tran- 
sition between CdGM's bound states is driven by the 
spectral flow although it is significantly modified by the 
impurities. This point of view provides us with some 
new knowledge about the dynamics of quasiparticles in 
the core of moving vortices. 

In this letter, in order to clarify the nature of the spec- 
tral flow, we first perform a finite size study of the en- 
ergy levels of quasiparticles in the presence of vortices. 
We consider a ring made from a strip of two dimensional 



(2D) BCS superconductor and concentrate on the case 
of single vortex in the system. The numerical evalua- 
tion of the eigenvalues of Bogoliubov-de Gennes equation 
clearly shows a spectral flow behavior. It is also shown 
that there are two kinds of spectral flow; one moving 
from the upper side to the lower side of the gap and the 
other moving in the opposite direction, which we call £+- 
and £_ -branch, respectively. 

Based on these facts, we argue that the transition of 
an electron from e_- to e + -branch, which occurs due to 
impurities, gives the microscopic description of the quasi- 
particle scattering, i.e. the dissipation in the vortex core. 
This process is studied from a quantum mechanical point 
of view employing Landau-Zener tunneling theory and it 
is clarified when the Bardeen-Stephen type description 
is valid. 

This letter is organized as follows: In Sec. || we de- 
scribe our model and basic scheme for the numerical cal- 
culation. Then the numerical result is introduced, based 
on which the nature of the spectral flow is discussed. In 
Sec. ||the effect of impurities is considered and the mech- 
anism of the dissipation in the vortex core is discussed. 
In Sec. | we discuss the relation between our results and 
the preceding studies. Some further problems are also 
mentioned. In Sec. || we summarize our results. 

§2. Bogoliubov-de Gennes equation 

In this letter, we consider a ring made from a strip 
of 2D BCS superconductor (Fig. 1), whose wave func- 
tions and energy spectrum of the quasiparticles are given 
by the Bogoliubov-de Gennes equation. The width and 
length of the strip are given by w and I, respectively. 
Magnetic flux $ is trapped in the cylinder to generate 
circulating supercurrent. We consider the case of single 
vortex located at R(t) = {R x , Ry) at time t. 

Bogoliubov-de Gennes equation is given by,E3 
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(2.1) 



* E-mail: hayashi@cmt.is.tohoku.ac.jp 



2 



Masahiko Hayashi 



where 



hf_ 

2m 



V 



he 



A -E F + V imp (r) (2.2) 



with m, e and c being the mass and charge (e < 0) of 
electrons, and the velocity of light, respectively; A = 
($/£) e x with e^ being unit vector in ^-direction; u n and 
v n are the wave functions of the quasiparticles and e n is 
their energy eigenvalue; Vi mp (r) is the impurity potential. 
We disregard the fluctuation of the gauge field in this 
letter. 





Fig. 1. (a) Picture of the model system which we consider in this 
letter. Magnetic flux <E> is introduced to generate circulating 
supercurrent in the ring. A vortex (indicated by a dot) crosses 
the ring and causes decay of the supercurrent. (b) Equivalent 
area in x-y plane. Lines at x = and x = I are identified. The 
position of the vortex is denoted by R(t). 
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Fig. 2. Infinite number of mirror vortices generated by the 
boundary conditions of the system (in the middle). 



First we solve Eq. (2.1) with Vi mp (r) = by intro- 
ducing an approximate solution for A(r,i), instead of 
solving it selfconsistently. The approximate solution is 
constructed in the following way. It is well known that 
the condition of no outgoing supercurrent at the bound- 
ary is satisfied by introducing mirror vortices. In the 
present case an infinite series of mirror vortices appears, 
as is seen from Fig. 2, which is considered to be two 
square vortex lattices with opposite flux overlapped with 
each other with a shift in y direction. Therefore we can 
employ the order parameter of vortex lattice to construct 
an approximate order parameter for the, present system. 
We use Eilenberger's order parametenia) for the vortex 
lattice and apply a correction to the. amplitude part in 
the sense of Clem's approximation,!!^ in order to make 



the amplitude constant except for the core region; 
<f>±{r) = e -^vTR y +wf ^ (- + T ji y + ^ 2wi) 



4>±l 

A <? 



(2.3) 



where #3(2, r) is the elliptic theta function; 'bar's denote 
lengths divided by I, e.g., x — x/l; c is a positive con- 
stant which determines the radius of the core; Ao is the 
equilibrium energy gap outside the vortex core; <f>*_ is the 
complex conjugate of </>_ . 

We assumed parabolic dispersion for the Fermions and 
took into account only the states whose energy E satis- 
fies Ep — hu>r> < E < Ep + hbjp, where Hu>d is a cut- 
off parameter. The matrix elements of the Hamiltonian 
(i.e., the matrix of the left hand side of Eq. (2.1)) were 



calculated using this basis and energy eigenvalues were 
obtained after diagonalization. 



In Fig. 

Ep = 10 3 x e , A = 10 2 x eo and hojp, — 300 x eo 



3, we show the result for the case of w — 0.5x1, 
> 3 x e A - - ln2 

■ 2 



where e = ( h A j2m) (2tt//) 2 . Here £ = hv F /A = 
n~ 1 l(Ep/ Ao)y/eo/Ep. Although the system is much 
smaller than the realistic one, our result shows a clear 
spectral flow behavior. 

Here we note a remarkable character of the spectral 
flow, namely the existence of many level crossings. This 
is the most essential process which leads to the dissipa- 
tion in the vortex core, as we see in the next section. This 
character is understood in the following way: If (u n , v n ) 
is an eigenfunction with eigenvalue e-_ K,«*J is also 
an eigenfunction with eigenvalue — e n L3 Therefore the 
energy spectrum is symmetric with respect to the Fermi 
energy. This means that if there are several spectral flow 
branches which are flowing down crossing the Fermi en- 
ergy, the same number of branches must be flowing up. 
These branches inevitably make crossing points. 

§3. Effect of impurity scattering 

In this section we study the effect of impurities 
^imp(r). For this purpose it is useful here to look at 
the dynamics of quasiparticles in superconductor-ncuyual 
metal-superconductor (S-N-S) Josephson junctionJl3't3' 
which shows a similar spectral flow behavior. a In S-N-S 
Josephson junctions, the bound states are formed in N- 
region due to Andreev scattering at S-N boundaries and, 
if we change the phase difference of two S"s, the eneijgy 
levels of the bound states also change continuously S$ 
We can also see level crossings similar with our system. 
When a finite bias voltage is applied to the junction, the 
phase difference increases in time as is known from the 
Josephson's acceleration relation and the spectral flow 
occurs. 

A precise analysis of this system, when N region is a 
point contact, was given by Averin and BardasJiZP They 
introduced a potential barrier in A-region and found 
that the scattering of quasiparticles by the barrier causes 
splitting of the crossing points of the spectral flow. The 
dominant scattering process is the backward scattering 
and the quasiparticles are scattered into the states with 
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Fig. 3. Energy spectrum of quasiparticles as a function of R y . We can see the energy gap opening at R y = and R y = w = 0.5 X I, 
and many energy level crossings in the middle. The crossing points are almost uniformly distributed in the range < R y < w except 
for the region which is approximately within the coherence length £ ~ 0.1 X / from the boundaries. The gap at R y = w is smaller 
than Aq = 10 2 X eo, because the average of the gap is suppressed by the effect of finite system size. The gap at R y = is suppressed 
further, because the supercurrent is larger at R y = than at R y = w. Some of the energy levels outside the gap are not plotted. At 
some points, splitting of the crossing points is seen, which is due to the approximation and not essential. 



opposite momentum. They also found that the dynamics 
of the scattering is describable employing Landau-Zener 
tunneling theory. 

From the analogy to the Josephson junctions, we can 
probably expect that the similar process occurs due to 
impurity scattering at each crossing points of Fig. 3. We 
can approximately estimate the magnitude of the split- 
ting as h/ri mp where r; mp is the life time of quasiparticles 
in the normal state due to impurity scattering. When 
the vortex is moving with a velocity vr, the scattering 
process can be described by the Landau-Zener tunneling 
theory. The time evolution of the system is shown in Fig. 
4 (a): If both of two crossing levels are vacant (or occu- 
pied) before crossing, the problem becomes rather triv- 
ial, since nothing special happens at the crossing point. 
However, if the lower level is occupied and the upper one 
is not, the final state changes depending on Vr. When vr 
is small, the particle follows the adiabatic change of the 
energy level, i.e., the broken curve in Fig. 4 (a). When 
vr is large, the particle has little time to interact with 
impurities and travels along the original levels without 
being scattered, i.e., dotted line. The probability of the 
latter case is given by 

where e+ — e_ is the relative velocity of two levels in 
energy space, which is estimated as follows: At low tem- 
peratures ksT < hivo, only the level crossings near the 
Fermi energy are relevant, which is depicted in Fig. 4 
(b). Firstly, the maximum of the spacing between two 
levels nearest to the Fermi energy is of the order of hujQ , 



which approximately correspond to the CdGM's bound 
states with angular momentum \x = 1/2 and —1/2. Sec- 
ondly, in our system the wave number of electrons in 
y-direction is quantized by tt/w. Therefore there exist 
kF /{tt/w) channels, each of which is responsible for one 
pair of spectral flow (e + and £_). If we assume that the 
crossing points near the Fermi energy are distributed uni- 
formly in < y < w when £ <C w, we can estimate the 
separation between crossing points to be n/kp, a half of 
the Fermi wave length. Although this is just an assump- 
tion at this stage, it seems to be a natural one, as one 
can see from Fig. 3. Therefore, noting R y — VRt, we 
obtain e + — e_ ~ UuiqUevr/'k . 
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Fig. 4. (a) Level crossing and Landau-Zener tunneling process, 
(b) The energy and length scale near the crossing point. 



We can see from Eq. (3.1) that when VRkpT? ujq <C 1, 



the quasiparticle follows the change of the energy level 
adiabatically, which means that the quasiparticle is scat- 
tered by the impurity. In contrast with this, when 
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VRkFT^pWa ^> 1, the quasiparticlc is not scattered but 
becomes a particle-hole pair after some level crossings, 
which costs approximately the energy of Ao. This energy 
cost may probably suppress the vortex motion by acting 
as a nondissipative force. Further studies are needed to 
clarify the detail of this process, however. 

It is also important to note that two different situ- 
ations arise from the relation between the level spac- 
ing hujQ and the magnitude of splitting h/Ti mp . When 
hujQ 3> S./Ti m p the situation is depicted in Fig. 5 (a). 
The energy levels oscillates with a period of a half of the 
Fermi wave length. In this periodic potential for 

the vortex may appear due to the change of the quasi- 
particle energy. If we simply sum up the energies of the 
quasiparticles in the bound states, the amplitude of the 
periodic potential becomes the order of Hujq , which is 
not negligible at lower temperatures ksT < Tlloq. On 
the contrary, when hu>o <C h/T lmp , the system recovers 
translational symmetry, as is seen from Fig. 4 (b). 

The dissipation of Bardeen-Stephen type is, therefore, 
observed when ujjfcF7f mp wo <C 1 and ujq <C 1/rjmp are 
satisfied at the same time. 
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Fig. 5. Schematic view of the energy levels of the bound states 
in the case of (a) u>o 3> 1/Ti mp and (b) u>o <C l/ri mp . 



§4. Discussion 

In the pioneering study, Caroli, de Gennes and Ma- 
tricon assumed that the vortex has rotational symme- 
try .HP Therefore the angular momentum was taken to be 
a good quantum number. This, however, is not the case 
with our system, in which the boundaries break the ro- 
tational symmetry and cause mixing of different angular 
momentum states. This is the origin of the spectral 
in this letter. Volovik, Makhlin and MisirpashaevJ 
by introducing the laboratory frame, obtained a similar 
spectral flow, which we consider is a realization of our 
finite size calculations in the bulk. The relation between 
two, however, is not clear at this stage. 

There are several works which have studied the dis- 
sipation in the vortex core from a quantum mechanical 
point of view, by addressing the transition of quasiparti- 
cles between diffemnt^dGM's bound states due to im- 
purity scattering S3 nl3'i3 The largest difference of our 
study from those is that the transition of quasiparticles 
is driven by the spectral flow which exists even with- 
out impurity scattering. In our treatment, however, the 
effect of impurity scattering is only introduced by Ti mp 
which is too simplified, especially in the case of small 



concentration of impurities. If we treat this more ex- 
actly, mesoscopic fluctuations may appear in the energy 
levels of Fig._5, reflecting the character of individual 
impurities Jilf E3 1 

In this letter we neglected the gauge fluctuation, 
namely the effect of finite penetration depth. Since the 
mirror vortices and their long-ranged phase modulation 
are important in our analysis, the effect of gauge fluctua- 
tion can be crucial. In order to clarify this point, we have 
to include the effect of gauge field into our Bogoliubov-de 
Gennes equation, which however needs further studies. 

In actual superconductors, three dimensionality is also 
important. In this case, the vortex become a line and 
each bound state in the vortex core obtains a continuous 
level corresponding to the motion along the line, which, 
of course, makes the problem much more complicated. 
However, since the dissipation is mostly due to the back- 
ward scattering of quasiparticles, which is expected to 
occur between the states with the same momentum in 
the direction of the line, we may believe that the phys- 
ical picture presented in this letter is still useful to un- 
derstand the basic mechanism of the dissipation. 

§5. Summary 

In this letter we studied the spectral flow and the ori- 
gin of dissipation in the core of superconducting vor- 
tices from a quantum mechanical point of view. We 
found that several characteristic regimes are defined by 
the material parameters and the vortex velocity. The 
dissipation of Bardeen-Stephen type occurs only when 
VRkFT? mp uJo <C 1 and loq <C 1/Ti mp are satisfied. In other 
regimes the effect of the spectral flow does not result in 
simple dissipation but can cause nondissipative forces as 
well. 
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